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CALCULATION  OF  HIGHER  APPROXIMATIONS  FOR  TWO-DIMENSIONAL 
COMPRESSIBLE  FLOW  BY  A  SIMPLIFIED  ITER  ATI  OF  PROCESS 
By  W.  H.  Braun  and  M.  M.  Klein 


SUMMARY 

The  iteration  equations  for  a  simplified  solution  of  the  nonlinear 
compressible -flow  equation  are  developed  and  applied  to  two  profiles, 
the  ellipse  and  the  Kaplan  section.  Limitations  imposed  on  the  solu¬ 
tions  by  simplifications  in  the  differential  equation  and  the  boundary 
conditions  are  discussed.  The  velocities  near  the  midchord,  the  criti¬ 
cal  and  potential  limit  free-stream  Mach  numbers,  and  the  extent  of 
isentropic  supersonic  regions  are  calculated  to  four  approximations 
for  the  Kaplan  section  and  to  six  approximations  for  the  ellipse  .  The 
development  of  the  iteration  equations  and  the  presentation  of  the 
results  are  made  in  conformity  with  the  Karman  transonic  similarity 
law  and  comparisons  are  made  with  other  solutions. 


INTRODUCTION 

Because  of  its  nonlinear  nature,  the  partial  differential  equa¬ 
tion  that  describes  the  flow  of  a  perfect  compressible  fluid  has  not 
yielded  completely  to  any  method  of  solution  so  far  advanced.  In 
fact,  in  such  methods  as  the  iteration  procedures  and  the  variational 
method  (reference  l) ,  it  is  not  even  attempted  to  find  a  closed-form, 
solution  or  the  general  term  of  a  series  expansion;  rather,  a  finite 
(and  usually  small)  number  of  terms  of  a  series  expansion  are  obtained. 
Unfortunately,  in  such  cases,  whereas  more  and  more  terms  are  required 
as  compressibility  effects  become  more  pronounced,  the  difficulty  of 
securing  the  higher  terms  becomes  inordinately  greater.  Although  the 
difficulties  of  nonlinearity  may  be  overcome  by  a  transformation  to 
the  holograph  plane  and  the  choice  of  a  convenient  gas  law  (refer¬ 
ence  2),  a  new  problem  arises  in  the  solution  of  any  boundary-value 
problem  because  of  the  distortion  of  the  boundary. 

The  solutions  by  iteration  are  begun  by  expanding  the  velocity 
potential  as  a  series  in  terms  of  some  flow  parameter;  it  is  then  pos¬ 
sible  to  determine  each  term  from  those  which  precede  it  by  solving  a 
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linear  equation.  Whenever  the  flow  is  over  a  thin  body  at  relatively 
high  speeds,  it  is  preferable  to  expand  the  potential  in  terms  of  some 
parameter  characteristic  of  the  body,  such  as  thickness  ratio,  camber 
coefficient,  or  angle  of  attack.  The  first  term  of  the  solution  is  a 
constant  representing  the  undisturbed  flow  and  the  second  is  the  well- 
known  linearized  or  small  perturbation  compressible  flow  (Prandtl- 
Glauert  solution) .  The  procedure  for  obtaining  successive  terms  is 
outlined  by  Ackeret  (reference  3)  and  extended  by  Kaplan  (reference  4), 
who  applies  it  to  a  bump  with  no  stagnation  points  and  to  a  circular - 
arc  profile.  The  higher-order  terms  were  shown  to  be  of  greater  con¬ 
sequence  as  the  Mach  number  increased. 

Recently  Perl  and  Klein  (reference  5)  applied  the  Prandtl- Ackeret 
procedure,  without  specifying  beforehand  the  parameter  of  expansion, 
to  flows  about  thin  bodies  at  transonic  speeds.  In  their  analysis, 
the  flow  equation  is  transformed  to  a  new  set  of  coordinates  in  which 
it  is  effectively  expanded  in  powers  of  (32  =  1  -  Mq2,  where  Mq  is 
the  free-stream  Mach  number.  In  the  transonic  range  ( P2  **  0) ,  the 
lowest  order  terms  in  p2  are  dominant ,  and  the  transonic  part  of  the 
flow  equation  is  easily  extracted. 

The  resulting  transonic  equation  does  not  of  itself  preclude-  size¬ 
able  perturbations  as  long  as  p2  is  chosen  sufficiently  close  to 
zero)  however,  the  Perl-Klein  process  uses  the  customary  small  pertur¬ 
bation  boundary  condition,  so  that  the  solution  of  the  complete  prob¬ 
lem  -  differential  equation  and  boundary  conditions  -  is  termed  the 
transonic  limiting  solution,  limiting  in  the  sense  of  Mach  number  near 
unity  and  geometric  parameter  near  zero .  This  limiting  nature  of  the 
solution  is  further  emphasized  by  the  form  of  the  expansion  parameter, 
which  is  found  to  be  simply  Karman!s  transonic  similarity  parameter 
(reference  6),  a  quantity  which  combines  the  three  definitive  character¬ 
istics  of  the  flow:  the  ratio  of  specific  heats  {that  is,  the  type 
of  gas),  the  free-stream  Mach  number,  and  the  geometric  parameter. 

Although  the  iteration  process  greatly  facilitates  the  solution 
of  a  nonlinear  problem  by  substituting  for  it  a  series  of  linear  prob¬ 
lems,  it,  in  turn,  raises  the  difficult  question  of  convergence  of  the 
series  solution.  The  present  analysis  was  made  at  the  NACA  Lewis  lab¬ 
oratory  to  extend  the  calculations  necessary  to  estimate  the  convergence 
of  the  iteration  process  of  reference  6  where  it  is  applied  to  two 
profiles,  the  ellipse  and  the  Kaplan  section.  In  particular,  it  will 
be  found  that  the  successive  approximations  suggest  convergence  not 
only  of  completely  subsonic  flows  but  also  of  some  isentropic  mixed 
flows . 

There  has  been  some  dispute  concerning  the  stability  of  the  second 
type  flow  (see  the  resume  of  these  discussions  by  Sears,  reference  7)) 
one  of  the  latest  contributions  is  a  stability  investigation  by  Kuo 
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(reference  8),  which,  indicates  that  decelerating  transonic  flows  are 
unstable  for  a  certain  type  disturbance.  The  result  of  the  instability 
is  the  formation  of  a  terminal  shock  in  the  supersonic  region,  with  a 
consequent  asymmetry.  Some  British  measurements  (references  9  and  10) , 
however,  indicate  that  small  supersonic • regions  not  terminated  by  shock 
waves  may  exist,  or  if  shocks  do  occur,  it  is  possible  that  they  may 
not  change  the  flow  greatly  from  the  isentropic  pattern,  provided  that 
the  free-stream  Mach  number  is  only  slightly  greater  than  critical. 

The  isentropic  mixed  flows  are  therefore  dealt  with  herein  under  the 
assumption  that  they  are  indicative  of  the  nature  of  some  real  flows. 

Another  unpleasant  feature  of  iteration  methods,  in  addition  to 
the  question  of  convergence,  is  the  large  amount  of  labor  involved.  In 
obtaining  the  higher  approximations  by  the  present  method,  it  is  con¬ 
venient  to  make  certain  simplifications  in  the  differential  equation 
so  that  the  labor  does  not  become  too  great.  The  justification  for 
such  changes  and  the  limitations  they  impose  on  the  solution  will  be 
discussed  as  they  arise. 


SMALL -PERTURBATION  TRANSONIC  FLOW 

In  reference  5  it  is  shown  that  the  partial  differential  equation 
governing  two-dimensional,  irrotational,  compressible  flow  may  be  writ¬ 
ten 

(  (3  (p^  +  <Pyy )  |i  Mq  (  2(PX  +  <PX  +  ^ 

=  [rM  (2<px  +  <px2)  +  r~  Mq4<1>/|<PXx  +  Mo2<Py2<Pyy  +  W  U  +  <PX)  Vxy 


(1) 


where 

P2  =  1  -  Mq2 

cp  perturbation  potential 
X  ratio  of  specific  heats  of  the  gas 
Mq  free-stream  Mach  number 

rM  =  MO2  d  +  ^Mq2) 


V 
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anri  the  free-stream  velocity  has  magnitude  unity  and  the  direction  of 
positive  x.  (A  complete  list  of  symbols  is  given  in  the  appendix.) 
If  the  flow  is  over  an  isolated  body  whose  surface  is  defined  by 

y  =  t  g(x) 

where  t  is  the  thickness  ratio,  the  boundary  conditions  on  the  per¬ 
turbation  potential  are 

<p  =  t(1  +  <px)  gx(x)  at  y  =  x  g(x) 


(Px  =  «Py  =  0 


at  x2  +  y2  =  ® 


For  purposes  of  solution,  it  is  desirable  to  write  the  second-order 
terms  of  equation  (l)  in  normal  form,  as  may  be  accomplished  by  the 
Prandtl-Glauert  transformation.  If  the  transformation  is  modified 
slightly,  however,  the  coefficient  of  a  nonlinear  term  of  the  same 
order  in  (32  also  becomes  independent  of  flow  parameters.  The  modi¬ 
fied  transformation  is 


x  =  x 


n  =  Py 


F(x,T])  =  -M  <p(x,y) 

P2 


and  it  transforms  the  problem  consisting  of  the  differential  equation  (l) 
and  the  boundary  conditions  (2)  into 


.  2  l  M02  1  9 

(1  -  2fx)fxx+  F^  =  P  \(r-l) rT-  F/F  +  JT  fxS 


2Mo2  \ 

+  r  ftfxt )/ 


'4  fa  *  r-f  p  F'12f-  +  2  P  W*) 

\  rM  rM  XM  / 


p6  /  lli  F  2AF  +  ^  F  S 

\  2  r  2  T]  p  2  T]  T]T| 

\  M  M  > 
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(1  +  <PV) 


at  Tj  =  PTg(x) 


Fx  =  =  0  at  x2  +  T)2/02  =  «  (6) 

p 

Figure  1  is  the  reciprocal  of  the  transformation  coefficient  rM/P  Of 
equation  (3) . 

A  transition  is  now  made  to  the  small -perturbation  transonic  prob¬ 
lem.  For  transonic  speeds,  the  parameter  p2  approaches  zero  so  that 
the  right  side  of  equation  (4)  may  be  neglected.  Because  of  the  non¬ 
linearity  of  the  right  side,  omission  of  these  terms  is  especially 
justified  if  the  perturbations  are  small.  If  the  body  considered  is 
thin  (l  >>t)  and  there  are  no  stagnation  points,  (px<  <1  and  the 
velocity  <px  may  be  neglected  in  boundary  condition  (5) .  For  flows  in 
which  stagnation  points  occur,  the  effect  of  these  stagnation  points 
is  assumed  negligible  in  the  region  of  interest  in  the  flow  field. 
Furthermore,  in  boundary  condition  (5),  because  the  product  p  is  a 
very  small  number  and  because  (as  will  be  shown)  the  solution  for 
equation  (7),  which  follows,  has  no  singularities  on  the  x-axis,  It 
is  assumed  that  the  condition  (that  the  velocity  normal  to  the  body  is 
zero)  may  be  made  on  the  x-axis  rather  than  on  the  body  itself.  This 
assumption  is  justified  in  reference  5.  Consequently,  the  small- 
perturbation  transonic  flow  problem  is  defined  by 

(l  2Fx)  FXX  +  FT]T]  ~  0  (7) 


ri  -z;  £>X 

p5 


T)  =  0 


F  =  F  =  0 

X  A  T] 


:2  +  T)2/P2  = 


For  values  of  Fx  <  i,  the  flow  equation  (?)  is  elliptic  and.  for 

U 


1. 

Fx  >  2  ^  kyPer^ol:i-c  •  In  figure  2  the  parabolic  value  of  the  per- 

turbation  velocity,  «PX  =  — —  (or  Fx  =  — ),  is  compared  with  the  excess 

2Fm  2 

of  the  critical  velocity  over  free-stream  velocity,  that  is,  the  pertur¬ 
bation  velocity  at  which  a  flow  becomes  locally  sonic,  as  calculated 
from  the  Bernoulli  equation.  For  such  values  of  M q  for  which  the  two 
curves  substantially  agree,  the  differential  equation  (7)  may  be  expected 
to  describe  mixed  flows  properly.  The  agreement  is  acceptable  in  the 
range  0.9  <  Mq  <  1.1,  which  is  approximately  the  range  of  accuracy  which 
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would  have  been  suggested  by  equabion  (4)  .  Ib  is  also  observed  ab  bhis 
poinb  bhab  if  bhe  error  in  bhe  boundary  condibion  (8)  is  bo  be  limibed 
bo  10  percenb  as  compared  wibh  bhe  exacb  boundary  condibion  (5) ,  ib  is 

necessary  bhab  <PX  <  Thus  an  examinabion  of  bobh  bhe  differenbial 

equabion  and  bhe  boundary  condibion  indicabes  bhab  flows  in  bhe .  bransonic 
domain  will  be  described  sufficienbly  accurabely  only  if  bhey  lie  in  bhe 

inberval  0< 

The  flows  bo  be  solved  herein  will  be  limibed  bo  bhose  aboub  bodies 
which  are  symmebric  aboub  bhe  x-  and  y-axes.  A  limibabion  bo  subsonic 
f ree-sbream  Mach  numbers  is  also  necessary  in  order  bhab  bhe  brans - 
formabion  (8)  does  nob-  become  imaginary.  The  resulbs  of  bhe  presenb 
mebhod  for  high  subsonic  speeds  can  be  applied,  however,  bo  slightly 
supersonic  free -sbr earn  flows  for  which  bhe  debached  shock  wave  is  far 
ahead  of  bhe  body  and  bhe  Mach  number  downsbream  of  bhe  shock  is  as 
much  below  uniby  as  bhe  Mach  number  upsbream  is  above  uniby.  In  bhis 
sense,  bhe  analysis  is  applicable  bo  bhe  enbire  range  of  bransonic 
velocibies . 


ITERATION  PROCESS 


For  convenience  in  describing  bhe  body ,  a  brans  formabion  is  made 
bo  ellipbic  qoordinabes,  one  of  which  is  cyclic  and  can  indicabe  posi- 
bion  around  bhe  circumf erence  of  bhe  body .  These  coordinabes  are  com¬ 
pared  wibh  carbesian  coordinabes  in  figure  3.  The  equabions  defining 
bhe  bransformabion  are 


x  =  cosh  s  cos  b 


rj  =  sinh  s  sin  b 


¥) 


=  cosh2  s  -  cos2  t  =  sinh2  s  +  sin2  t 


(10) 

(11) 

(12) 


where  J 


(s) 1 


is  the  Jacohian  or  functional  determinant  of  the 


x,y-coordinates  with  respect  to  s  and  t.  The  reciprocal  of  the 
Jacohian  of  the  transformation  is  plotted  in  figure  4  as  a  function 
of  t  for  several  values  of  s . 

The  first  derivatives  of  the  potential  in  the  new  system  are 
Fx  =  (sinh  s  cos  t  Fg  -  cosh  s  sin  t  Ft)J_1 
F  =  (cosh  s  sin  t  Fg  +  sinh  s  cos  t  Ft)j_1 


(13) 

(14) 
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which  become  at  the  axis  of  the  body 


F  =  -  — - —  F+ 
x  sin  t  "t 


Fn  sin  t 


(15) 

(16) 


Substitution  of  these  into  the  boundary  condition  (8)  yields 


tF 


F  = 
r  s 


M 

3 


s  =  0 


(17) 


Boundary  condition  (9)  becomes 


* 


4 


V 


Fs  =  Ft  =  °  s  =  •  (18) 

An  exact  transformation  of  the  flow  equation  (7)  would  result  in 
the  inverse  Jacobian  appearing  in  the  nonlinear  term,  and  it  would  be 
necessary  to  expand  it  in  a  Fourier  series  before  employing  a  method 
of  solution  involving  the  separation  of  variables.  The  process  used 
herein  would  then  require  the  solution  of  an  infinite  number  of  ordinary 
differential  equations.  In  order  to  simplify  the  analysis  and  to  reduce 
the  computational  work  that  would  ordinarily  accompany  such  an  expan¬ 
sion,  only  its  first  or  constant  term  is  retained.  That  is,  in  the 
transformation  of  equation  (7)  to  elliptic  coordinates,  the  Jacobian 
is  set  equal  to  1,  which  is  the  first  term  in  its  expansion  about  the 
point  s  =  0,  t  =  jt/2.  Of  course,  this  simplification  will  restrict 
the  region  in  which  the  equation  is  valid.  A  further  restriction  on 
the  size  of  the  region  of  validity  is  made  by  restricting  the  calcula¬ 
tions  to  velocities  on  or  near  the  surface  of  the  body;  thus  s  *  0  in 
the  differential  equation. 

If  the  relation 


■^xx  +  ^tjt]  J  (Fss  +  Ftt) 

and  equations  (13)  and  (14)  for  the  derivatives  of  Fx  and  F~  are 
used,  equation  (7)  with  the  approximations 

J  *  1 


and 


sinh  s  *  0 
cosh  s  *  1 
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may  be  -written 

Fss  +  Ftt  =  -  2  ain2t  Ft  ~  (sin  t  Ft) 

=  -  2  sin^t  F^F.^  -  2  sin^t  cos  t  F^-F^.  (19) 

This  equation  is  proposed  for  use  as  a  description  of  the  flow  in  a 
region  about  the  origin  whose  boundary  is  set  by  the  criteria 

j-1  _  pi  <  0.10  and  s  <  0.1.  Such  a  region  in  the  xT}-plane  is  shown 
in  figure  5,  where  for  comparison  a  10-psrcent -thick  (in  the  xy-plane) 
Kaplan  section  is  included  as  a  typical  profile.  The  length  of  the 
rectangle  is  seen  to  be  about  0.3  chord. 

With  the  change  in  the  region  of  validity  of  the  differential 
equation,  there  must  be  a  change  in  the  boundary  conditions  to  corres¬ 
pond  to  this  region.  The  surface  boundary  condition  (7)  may  be 
retained  as  is;  therefore  only  the  outer  boundaries  (dashed  rectangle 
in  fig.  5)  of  the  region  lack  boundary  conditions.  As  a  substitute 
for  conditions  on  these  boundaries,  the  boundary  condition  (18)  at 
infinity  is  taken  as  a  guide  to  postulate  that,  in  the  restricted 
region,  both  Fs  and  F-j-  must  decrease  with  increasing  s.  The  fav¬ 
orable  comparison  with  more  complete  solutions  shows  that  this  con¬ 
dition  prescribing  monotonic  character  of  the  potential  derivatives  is 
sufficient.  These  comparisons  are  made  in  subsequent  sections. 

The  method  of  solution  of  the  problem,  as  formulated  by  the  dif¬ 
ferential  equation  (19),  the  boundary  condition  (17),  and  the  condition 
of  monotoneity  in  s,  is  to  write  the  transformed  potential  F  as  an 
infinite  series 

12  3 

F=F+F+F+.  .  .  (20) 

where  it  is  assumed  that  the  succeeding  terms  become  progressively 
smaller  at  such  a  rate  as  to  insure  convergence  of  F  and  its  deriva¬ 
tives.  Definition  (20)  is  placed  in  the  flow  equation  (19),  which  is 
then  arranged  as  follows : 
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1  1 
Fss  +  Ftt 


+  Fss  +  Ftt 


-  2  sin3t  F^F.^  -  2  sin2t  cos  t  F^F^ 


+Fss  +  Ftt 


2  112  12 
2  sin3t  (FtFtt  +  F^F.^)  -  2  sin2t  cos  t  2  F^F. 


+Fss  +  Ftt 


n-1  v  n-v 


n-1  v  n-v 


-  2  sin3t  2  FtFtt  -  2  sin2t  cos  t  2  F^F. 


This  arrangement  indicates  that  the  equation  should  he  divided  for  pur¬ 
poses  of  iteration  into  a  sequence  of  equations  each  corresponding  to 
one  line  of  equation  (2l).  In  this  way  all  terms  which  have  the  same 
total  superscript  value  are  grouped  together;  this  procedure  will  later 
be  seen  equivalent  to  separating  out  terms  of  the  same  order  in  the 
Karman  similarity  parameter . 

The  boundary  condition  on  the  body  is  also  divided  in  a  particular 

1 

way  among  the  terms  in  the  series  solution.  The  first  term  F  makes 
the  total  contribution: 


Fs  = 


1 M 

.3  St 


s  =  0 


whereas  the  boundary  conditions  for  all  other  terms  are  homogeneous: 


n 

Fg=0  s  =  0  n  >  2  (23) 

n  n 

The  condition  of  monotoneity  requires  that  all  F  and  F.  are  monotone 

S  b 

decreasing  in  s. 

1 

FIRST  APPROXIMATION,  F 

The  first  approximation  differs  from  subsequent  approximations  in 
that  it  requires  the  solution  of  the  homogeneous  Laplace  equation  (first 
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line  of  equation  (2l))  with  a  nonhomogeneous  boundary  condition  (22). 

Higher  approximations  involve  the  solution  of  nonhomogeneous,  or  Poisson, 
equations  with  homogeneous  boundary  conditions.  For  this  reason  and 
to  obtain  the  necessary  start  in  the  iteration  procedure,  the  first 
approximation  is  treated  separately . 

The  solution  of  this  first  iteration,  which  takes  account  of  the 
monotoneity  of  the  derivatives  and  of  t  as  the  cyclic  coordinate,  is 

1  00 

F  =  2  e“kS  (ok  cos  kt  +  sin  kt)  (24) 

k=l 

Before  this  form  is  introduced  into  the  boundary  condition  (22),  the 
function  g  (0,t)  is  expanded  into  a  Fourier  series  and  written 

g  (0,t)  =  2  sin  jt  (25) 

j 

In  conformance  with  the  restriction  to  bodies  which  are  symmetric  about 
the  y-axis  and  about  which  there  is  no  circulation,  no  cosine  terms 
appear  in  the  expansion.  The  index  j  is  limited  to  odd  values 
because  of  the  symmetry  of  g  about  t  =  .2.  The  introduction  of  expres¬ 
sions  (24)  and  (25)  into  the  boundary  condition  on  the  body  (22)  leads 
to 

oo  "iTw  00 

-  2  k(ak  cos  kt  +  sin  kt)  =  -  —  2  j8j  cos  jt  (26) 

k=l  P  0=1 

j  odd  only 

It  follows  that 

F  =  2  Sk  e_]£S  cos  kt  (27) 

P5  k=l 

k  odd  only 


The  nature  of  the  dependence  of  F  on  ,  a  quantity  combining 

P6 

the  entire  dependence  of  the  first  approximation  on  the  parameters  of 
the  flow  (reference  6),  has  the  consequences  immediately  observed  from 
equation  (21)  that 
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(28) 


The  quantity 


j 


is  recognized  as  the  similarity  parameter  of  reference  5,  and  reduces 
to  Karman!s  similarity  parameter  (reference  6)  as  Mq  -►  1.  Hence ,  as 
previously  claimed,  equation  (2l)  is  the  ordering  in  powers  of  K  with 
respect  to  this  parameter.  The  similarity  parameter  is  plotted  as  a 
function  of  M q  for  various  %  in  figure  6,  which  is  reproduced  from 
reference  5. 


n  n 


n 


Writing  F  s  f K  gives  the  perturbation  potential,  by  the  trans¬ 

formation  (3),  as 


T  l1  2  3  2  \ 

<p  =  ^  \f  +  f K  +  flv  +  .  .  .) 


SX ft 

pL\  k=i 


■ks 


) 


.3_2 


e  cos  kt  I  +  f  K  +  f  K  + 


(29) 


The  first  term  inside  the  parenthesis  represents  an  incompressible  flow. 
Hence,  for  very  small  values  of  T ,  for  which  K  is  also  small,  the 
transonic  potential  (29)  reduces  to  the  Prandtl-Glauert  rule. 


HIGHER  APPROXIMATIONS 
n 

Inasmuch  as  all  the  F  beyond  the  first  satisfy  the  same  general 
form  of  differential  equation  and  have  the  same  boundary  condition,  it 


V 
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is  possible  to  treat  the  general  form  of  F  rather  than  each  approxi¬ 
mation  separately.  In  the  following  discussion,  an  induction  argument 
will  be  used  to  show  that  every  term  in  the  series  for  the  transformed 
potential  has  the  form 


X  ( |a+3 )  X  ( |i+3  ) 

X  ,  -3  -3 

F  =  n  2  2 

r=l  P=1 


XX  X  p  -ps  , 

(L^  +  MppS  +  N^s*  +  •  •  •  )e  *  cos  rt 

(30) 

r  odd  only 


XX  X 

where  Lrp,  Mrp,  and  Nrp  are  constants  and  (a  is  the  number  of  terms 

in  the  expansion  (equation. (25))  of  q(t),  the  function  characterizing 
the  shape  of  the  body.  The  main  part  of  the  argument  is  to  show  that 
if  equation  (30)  is  true  for  all  X  <  n,  then  it  is  also  true  for 

X  =  n. 


From  equation  (2l),  the  differential  equation  satisfied  hy  any 
n 

term  F  (n  >  l)  is 


n  n  n-l  v  n-v  n_l  v  n-v 

F  +  F. .  =  -  2  sin3t  2  FtFtt  -  2  sin2t  cos  t  2  FtFt  (3l) 

SS  tt  M  _l  V  =1 


* 


V 


v  n-v 

Use  of  the  postulated  form  (30)  for  F  and  F  leads  to  the  expan¬ 
sion  for  (3l) , 

9 


* 
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v 


* 


n 

ss  +  ^tt 


n-1 


v(^+3)  (n-v)  (p.+3)  v  ( |i+3 )  (n-v)(ki+3) 


.E  2  2 

8  t  i1 

v  =1  m=l 


-3 

2 

P=1 


-3 

2 

k=l 


-3 

2 

r=l 


n-v  n-v 


,n- v 


(Lkm  +  sJW  +  s  wkm  +  •  •  •)  (Lrp  +  sM^  +  s%p  +  .  .  .) 

—  ( ~p4-Tn )  s  I  o  o 

e  l(-3kr  -kr)  cos  (k+r+l)  t  +  (-3kr  +kr)  cos  (k-r+l)t  + 

(3kr^+kr)  cos  (k-r-l)t  +  (3kr^-kr)  cos  (k+r-l)t  + 


2 

(kr+kr  )  cos  (k+r+3)t  +  (kr-kr2)  cos  (k+r~3)t  + 
(-kr+kr^)  cos  (k-r+3)t  +  (-kr-kr2)  cos  (k-r-3)tj 


(32) 


k,  r  odd  only 


The  right  side  of  the  foregoing  equation  is  the  sum  of  products  involv- 

V  V  V 

ing  the  constants  ^km^  •  •  •>  powers  of  s,  exponentials  in 

s,  and  cosines  of  angles  of  multiple  t.  Accordingly,  after  all  five 
summations  have  heen  accomplished  and  regroupings  made  according  to 
powers  of  s,  exponentials  in  s,  and  multiples  of  t,  the  equation 
will  be 


n  n 


n(p+3)  n(p+3) 

-3  -6  n 


n  n 


Fss  +  Ftt  =  Rn  3,  +  sBh.1  +  s  ChJ  -  •  •  •)  *  =°»  ht 


h=l  j=2 


(33) 


h  odd  only 


If  a  solution 


n(n+3) 

n  -3  n 

F  =  K11  2  qh  (s)  cos  ht  (34) 

h 


V 


h  odd  only 
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is  assumed,  the  coefficient  of  each  harmonic  term  in  (33)  must  vanish 
Thus 


2  %.  (s)  "  h  %.  (s) 


n(p+3) 

-6  n  n  2n 

2  (Ahj  +  +  s  ChJ  + 

0=1 


.)  e"J; 


h  odd  only 

The  particular  integral  of  this  equation  is 

n  , 

A  ~hS 

n  n  e-JS  Aj^se 


.2 

_  2 

3 

-h 

p 

-js 

se 

.2 

,2 

J 

-h 

“js  5  -hs 


j  c  v£L  sz  +  JJjl.  +  £jW_  .  +  si  +  -S- 

hJj2-h2L  Ah2  m  v»  ^  ^ 


Here  the  notation  2  refers  to  a  summation  over  all  values  of  j 

j/h  n  -hs  11 

except  j  =  h.  The  complementary  function  is  G^e  ,  where  may 

he  found  from  the  boundary  condition  on  the  body  (23): 


^♦(S)  -o 

\ds  /s=o 


n  >  2 


Then 


S/77 
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where . 


n 

^bh  = 

« 


n 


n  _ i 

^  h  Us  J  s=0 


n 


e  +  Pi 


-he 


n(n+3) 

-6  n  n  pn  _n 

Z  (Ly  +  sM^  +  sXj  +  s5°hj  + 
j=l 


.  .)  e-Js 


(38) 


n 

Lh. 


=  2 


j/h  ^2-h2 


n  n 


■] 


Afr  Hi  -if) 


(n  n  n  \ 

lAhh  -^2  +  Bhh-^3+chh-^4+-  •  J 

\  2h  4h  4h  / 


+21i  +  4  j 


h2)2  (j2-h2J 


+  «  • 


n  ,  /n  4  •  n  \ 

Mhj  =  . j  .2  2  Whj  +  .2  ^2  Chj  +  ’  ’  ’  J 

jfh  j  -h  \  j  -h  / 

n  /nlnini  \ 

Mhh  =  "  \Ahh  2h  +  ^h  — 2  +  CKh  2  +  *  '  ’] 
\  4h  4h  1 


n 

K  . 
hj 


=  Z 


1  n 

C,  .  +  .  •  • 


j/hj2.h2 


N- 


hh 


n  i  n 


“  IBhh  TiT  +  Chh  ~  + 


4h 


) 


bh 


n  i 

-  c. .  +  . 

hh  6h 


(39) 
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(The  maximum  value  of  h  for  any  n  in  equations  (38)  and  (39)  Is  ^ 
n(|i+3)-3.)  Finally,  if  equation  (38)  is  placed  in  the  assumed  form  (34), 

the  n^*1  approximation  becomes 


F  =  K11 


n(p+3) 

-3 

2 

h=l 


n(p+3) 

-3  /n  n  on  3n 
2  Uhj  +  sM^  +  s  Nh<j  +  s  Otj  + 

0=1  x 


cos  ht 


(40) 


h  odd  only 


which,  is  exactly  the  form  postulated  in  equation  (30) . 

v 

It  has  "been  shown  that  if  F  has  the  form  (30)  for  all  v  <  n;  it 
will  also  have  that  form  for  v  =  nj  there  remains  only  to  test  this 
form  on  a  particular  approximation  in  order  to  complete  the  induction 
proof.  If  the  second  approximation  is  taken  to  correspond  to  the 
nth  approximation  of  the  proposition,  a  comparison  of  all  previous 
approximations  (equation  (27) ,  in  this  case)  indicates  that  the  n 
(second)  and,  hence,  all  approximations  subsequent  to  it  must  obey 
the  form  (30).  The  fact  that  the  method  of  solution  for  the  first 
approximations  differs  from  that  of  subsequent  approximations  does 
not  interfere  with  the  argument. 


The  actual  mechanics  of  the  iteration  process  have  been  derived  in 
the  foregoing  development  and  are  outlined  by  equations  (32),  (33), 

(39),  and  (40)  .  For  example,  after  the  first  approximation  has  been 
obtained,  the  differential  equation  for  the  second  may  be  found  by 
setting  n  =  2  in  equation  (32)  and  letting  k,  r,  m,  and  p  assume 

1 

values  in  accordance  with  the  form  of  F .  After  all  the  coefficients 
of  each  harmonic  on  the  right  side  have  been  gathered,  the  differential 
equation  will  have  the  form  (33).  The  application  of  the  formula  (39) 

2  2 

yields  the  constants  of  the  solution  (40)  in  terms  of  the  A^,  -^hj> 

112  2 

.  .  .  of  the  differential  equation.  The  Lhj,  •  •  •  >  ^hj; 

3  3 

.  .  .  are  then  known  and  are  sufficient  to  determine  the  L-^j,  Mj^jj  .  .  . 
if  the  foregoing  process  is  repeated. 


MAXIMUM  VELOCITY  INCREMENT  AND  CONVERGENCE 

For  flow  over  a  profile  symmetric  about  the  y-axis,  the  maximum 
velocity  increment  occurs  at  the  midchord.  At  the  critical  free-stream 
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4 


Mach  number,  the  flow  will  he  sonic  at  the  midchord,  and  for  higher 
free -stream  Mach  numbers,  a  local  supersonic  region  will  develop. 

Connected  with  the  formation  of  this  region  is>  the  question  of 
convergence  of  the  series  obtained  from  the  iteration  process  to  repre¬ 
sent  the  maximum  velocity  increment.  This  question  of  convergence  is 
especially  important  because  the  maximum  velocity  increment  is  expected 
to  be  the  first  velocity  component  at  any  point  in  the  flow  to  diverge, 
and  this  divergence  signals  the  failure  of  the  iteration  process  to 
serve  as  a  method  of  solution.  If  several  terms  of  a  series  are  known 
and  subsequent  terms  are  assumed  to  decrease  in  the  same  manner ,  an 
estimate  of  the  domain  of  convergence  may  be  made  as  follows : 

The  limiting  (t  -+  0,  x  0)  perturbation  velocity  at  the  midchord 
is  a  sum  of  terms  (equations  (3),  (15),  and  (30)) 


n 

lim  q>v 
x  -  0 

y  -*  0 

h-l 

=  |Kn-!  Z  ZLhj  h(-l)  2 
h  j 

s  |  K31-1  p  (41) 


or 


lim  <Px  =  1  Z  K11"1  P 
x  -►  0  H  n 

y  -*•  0 


(42) 


The 


are  seen  by  definition  to  depend  only  on  the 


hj 


which  intro¬ 


duce  the  effect  of  the  shape  of  the  body.  The  series  (42)  will  converge 
if  the  limit  of  the  ratio  of  successive  terms  is  less  than  1  for  all  n: 


n+1 

,  .  K11 

lxm  - ^  <  1 

n  -*•  •  x  J1-1  n 
D 


V 
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■» 


Kj  =  g.l.b.  — ^ —  n->®  (43) 

P 

where  Kz  is  the  maximum  value  of  the  similarity  parameter  for  which 
convergence  of  the  iteration  process  will  occur. 

The  so-called  potential  limit  Mach  number,  or  free-stream  Mach 
number  for  which  the  solution  diverges  for  any  T ,  is  defined  by 


Kz 


or 


(i  - 


K. 


(44) 


« 


» 


SURFACE  VELOCITIES 

For  small  disturbances  such  that  the  squares  of  the  perturbation 
velocities  may  be  neglected  compared  with  unity,  the  transformed  limit¬ 
ing  perturbation  velocity  on  the  ellipse  is 


F*  -  -  Hirt  Ft  (0’t) 


n 


sin  t 


Kn  2 


n 


2  h  sin  ht 

h  j 


or,  in  the  x,y-plane, 

<PX  (x'°)  =  2  K11"1  2  2  Lhj  h  sin  ht  (45) 

n  h  j 

Equation  (45)  will  hold  only  on  that  portion  of  the  surface  along  which 
•J  =  1  (figs.  5  and  6). 


<sr 


2275 


2275 


MCA  TN  2511 


19 


The  local  Mach  number  at  the  midchord  is  given  for  the  small  per¬ 
turbation  case  by  equation  (27)  of  reference  11: 

M2  -  1  =  -  32  +  2TM  <px  (0,0)  (46) 

In  a  previous  discussion  concerning  figure  2  it  was  seen  that  the 
hyperbolic  region  of  the  limiting  equation  (7)  represents  very  nearly 
the  region  of  supersonic  flows  for  a  perfect  fluid.  Therefore  the 
parabolic  curve 

2 

FX  =  |  or  **  =  ff~  (47) 

c  M 

can  be  used  to  outline  the  supersonic  region  or  to  find  significant 
points  on  it  in  the  physical  plane.  Whether  or  not  an  iteration  solu¬ 
tion  can  correctly  describe  a  supersonic  flow  may  be  questioned,  for 
as  performed  here  the  iteration  procedure  is  essentially  the  solution 
of  a  set  of  elliptic  equations .  The  applications  which  follow,  although 
not  answering  this  question  directly,  do  suggest  convergence  of  the 
solution  in  the  hyperbolic  region. 

For  any  thickness  ratio  t  the  critical  Mach  number,  or  lowest 
free-stream  Mach  number  at  which  some  point  in  the  flow  field  becomes 
sonic,  is  given  by 


rMT 


3  =Kc 


or 


(l  -  ¥  %2) 

(1  -  m//2 


(48) 


APPLICATION  TO  ELLIPTIC  CYLINDER 

The  expression  that  describes  the  upper  surface  of  an  ellipse 
whose  major  axis  coincides  with  the  direction  of  flow  is 


y  =  T  g(x)  =  x/\/l  -  x2 
=  t  sin  t  (s  -►  0) 


(49) 
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It  follows  from  equation  (27)  that 

F  =  Ke“S  cos  t  (50) 

Subsequent  approximations  through  the  sixth  have  been  calculated 
by  means  of  the  iteration  process  previously  described.  The  coefficients 

V  and  %  ^  the  V  %'  %  calculated  ^ 

means  of  equations  (33)  and  (39)  are  presented  in  table  I  in  matrix 
arrangement.  From  equation  (4l)* 


<P  (0,0)  =  \  (1  +  0.35237  K  +  0.36986  K2  +  0.51836  K3  + 
x  J3 

0.83711  K4  +  1.42014  K5  +  .  .  .  )  (5l) 

For  comparison*  the  maximum  disturbance  velocity  as  determined  by 
the  first  three  approximations  of  the  iteration  solution  herein  des¬ 
cribed  is  plotted  (fig.  7)  with  the  exact  solution  of  Hantzsche  (ref¬ 
erence  12).  Hantzsche *s  iteration  process  is  based  on  an  expansion  of 
the  stream  function  in  powers  of  the  thickness  ratio.  The  agreement 
between  the  two  methods  appears  close  enough  to  justify  the  calcula¬ 
tion  of  higher  approximations  for  obtaining  a  better  estimate  of  the 
potential  limit  Mach  number. 

An  estimate  is  made  in  figure  8  of  the  largest  value  of  the 

similarity  parameter  for  which  the  series  for  the  maximum  velocity 
increment  will  converge.  The  broken  line  joins  successive  estimates 
of  K i  obtained  by  comparing  terms  by  means  of  a  ratio  test  described 
in  the  section  "MAXIMUM  VELOCITY  INCREMENT  AND  CONVERGENCE."  An 
examination  of  the  slopes  of  the  segments  of  the  line  indicates  that 
the  limiting  ratio  of  these  slopes,  for  large  n,  is  about  l/2.  If 
it  is  assumed  that  these  slopes  continue  to  decrease  in  a  ratio  l/2, 
the  greatest  lower  bound  of  the  broken  line  curve  can  be  calculated 
by  a  geometric  progression.  The  estimated  greatest  lower  bound  of  the 
broken  line  is  then  about  =  0.56. 

That  the  solution  converges  even  when  some  of  the  flow  is  super¬ 
sonic  (hyperbolic)  is  indicated  in  figure  9,  which  shows  the  contribu¬ 
tion  of  each  approximation  to  the  maximum  velocity  increment  for  vary¬ 
ing  K.  Supersonic  flows  lie  to  the  right  of  the  parabolic  curve  and 
the  intersection  of  this  curve  and  the  sixth  approximation  marks  the 
value  of  the  parameter  K  at  which  the  transition  from  totally  sub¬ 
sonic  flows  to  partly  supersonic  flows  occurs.  This  transition 
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point,  Kc  =  0.595,  is  also  very  evident  in  figure  10  as  the  intersection 
of  the  curves,  for  several  values  of  M q,  of  the  local  Mach  number  at 
midchord  (equations  (42)  and  (46)). 

The  surface  velocity  on  the  ellipse  for  three  values  of  K  is 
presented  in  figure  11.  The  parabolic  line,  equation  (47),  is  included 
in  each  case  for  comparison.  Thus  in  figure  11(a)  all  velocities  are 
subsonic,  whereas  figure  ll(b)  (K  =  0.4)  shows  a  partly  supersonic 
flow.  This  is  in  agreement  with  the  previous  calculation  showing 
transition  from  subsonic  flow  to  supersonic  flow  at  Kc  =  0.395.  The 
longitudinal  extent  of  the  symmetric  supersonic  region  is  seen  to  be 
about  l/lO  chord,  based  on  six  approximations.  The  last  of  this  series 
of  figures  (fig.  11(c))  is  for  a  value  of  K  near  .  Although  the 
series  for  the  velocity  increment  may  be  converging,  it  is  evident 
that  near  the  midchord  the  uncalculated  approximations  are  not  neglig¬ 
ible  for  this  value  of  K.  The  width  of  the  supersonic  region  is 
greater  than  the  extent  of  the  region  of  accuracy  of  the  calculations. 

Curves  of  critical  Mach  number  and  potential  limit  Mach  number 
against  thickness  ratio  have  been  plotted  in  figure  12  using  equa¬ 
tions  (48)  and  (44),  respectively. 


APPLICATION  TO  KAPLAN  BUMP 


As  a  second  example,  the  iteration  process  is  applied  to  the  cal¬ 
culation  of  noncirculatory  flow  past  a  Kaplan  section;  this  problem 
has  been  treated  in  an  exact  manner  to  three  approximations  in  refer¬ 
ence  4.  The  Kaplan  section  has  no  stagnation  points  (fig.  5)  so  that 
nowhere  is  the  boundary  condition  singular.  Its  parametric  repre¬ 
sentation  is  (reference  4), 


x  =  cos  t 

y  =  j  (3  sin  t  -  sin  3t) 
Application  of  equation  (27)  yields 


1 

F  = 


cos  t  - 


-3s 

e 


cos  3t) 


(52) 


(53) 


The  second,  third,  and  fourth  approximations  have  been  calculated 
according  to  the  procedure  given  in  the  section  "HIGHER  APPROXIMATIONS" 

n  n 

and  are  recorded  in  table  II  as  arrays  of  the  constants  Lua,  Mu  4,  and 
n  n  n  n 


N^j .  The  constants  B^j,  and  C^j 

terms  of  the  differential  equations  (33). 


represent  the  nonhomo gene ous 
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The  chord-wise  perturbation-velocity  distribution  at  the  surface  is 
plotted  in  figure  13  for  comparison  with  the  result  of  Kaplan  (refer¬ 
ence  4)  for  t  =  0.10  and  Mq  =  0.75.  The  agreement  is  probably  suf¬ 
ficient  from  a  practical  standpoint  to  justify  the  labor-saving  limit¬ 
ing  processes  that  have  been  made  in  the  iteration  process. 

An  estimate  of  the  limiting  value  of  the  parameter  K  for  which 
the  series  for  the  maximum  velocity  increment  will  converge  is  made 
in  figure  14,  where  the  ratios  of  successive  coefficients  in  the  series 
for  the  maximum  velocity  increment  are  plotted  for  the  first  four  terms . 
The  series  for  (<Px)max  'to  four  approximations  is  calculated  from 
equation  (4l)  to  be 

(<p  )  =  +  0.81762  K  +  1.15715  K2  +  2.11564  K3  +  .  .  .)  (54) 

v  x/max  B  J 


An  analysis  similar  to  that  employed  on  the  corresponding  series  for 
the  elliptic  cylinder  indicates  that  the  series  (56)  ‘will  converge  for 
values  of  K  <  0.387. 

The  value  =  0.387  is  considerably  above  the  value  of  K  for 
•which  the  flow  becomes  sonic  at  the  midchord,  as  given  in  figure  15  by 
the  intersection  of  the  fourth  approximation  curve  and  the  parabolic 
curve.  The  converging  nature  of  the  successive  approximations  indi¬ 
cates  that  this  value  of  Kc  =  0.270  will  be  lowered  only  slightly 
by  increased  accuracy  from  higher-order  terms. 

The  variation  of  the  local  Mach  number  corresponding  to  the  maxi¬ 
mum  velocity  increment  has  been  calculated  by  means  of  equation  (46) 
and  is  plotted  against  the  similarity  parameter  for  several  values  of 
free-stream  Mach  number  in  figure  16.  The  calculations  are  based  on 
four  approximations .  The  velocity  increment  on  the  surface  of  the  bump 
is  plotted  as  a  function  of  distance  along  the  chord  in  figure  17  .  Of 
these,  figure  17(a)  shows  the  only  completely  subsonic  flow.  In  fig¬ 
ure  17(b),  Kc  <  K  =  0.3  <  Kj.  Hence  the  flow  is  everywhere  convergent, 
but  near  the  midchord  there  is  a  small  supersonic  region  represented  by 
those  values  of  (Fx)max  above  the  parabolic  curve.  For  figure  17(c), 
the  similarity  parameter  is  chosen  slightly  larger  than  the  estimated 
maximum  for  convergence,  as  obtained  from  equation  (43)  and  fig¬ 

ure  14;  and,  although  successive  approximations  still  decrease  monoton- 
ically,  the  curves  of  figure  17(c)  should  not  be  expected  ultimately 
to  converge. 

Figure  18  shows  the  critical  Mach  number  Mq  and  the  potential 
limit  Mach  number  Mj  as  functions  of  thickness  ratio.  For  compari¬ 
son,  Mc  and  M^  as  computed  by  Kaplan  (reference  4)  are  included. 

The  present  calculations  indicate  that  isentropic  mixed  flow  is 
restricted  to  a  considerably  smaller  range  of  Mach  number. 
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CONCLUDING  REMARKS 

An  iteration  procedure  similar  to  that  employed  in  the  Ackeret- 
Prandtl  type  procedure  has  been  utilized  in  the  calculation  of  the 
symmetric -type  transonic  potential  flow  for  an  elliptic  cylinder  and 
a  Kaplan  section.  In  order  to  obtain  an  estimate  of  the  region  of  con¬ 
vergence  of  this  type  of  solution,  the  work  has  been  carried  through 
six  approximations  for  the  elliptic  cylinder  and  through  four  approxi¬ 
mations  for  the  Kaplan  section. 

The  results  indicate  that  the  iteration  solution  converges  for  a 
range  of  the  similarity  parameter  for  which  supersonic  regions  exist 
in  the  flow  field.  For  these  hyperbolic  regions,  the  nonlinear  term 
of  the  potential  equation  is  not>  strictly  speaking,  small  compared 
with  either  term  of  the  Laplacian.  The  iteration  procedure  presented 
herein  may  therefore  be  applicable  to  a  greater  range  of  disturbance 
velocities  than  would  be  indicated  by  its  method  of  formulation,  in 
which  the  nonlinear  term  is  considered  a  perturbation  on  the  Laplacian. 

The  results  have  been  presented  in  conformity  with  the  transonic 
similarity  law.  The  essential  characteristics  of  a  flow  may  then  be 
presented  in  a  simple  and  direct  manner  by  classifying  the  flow  with 
respect  to  significant  values  of  the  similarity  parameter,  for  example, 
the  critical  and  potential  limit  values . 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  June  26,  1951 
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APPENDIX  -  SYMBOLS 


Tiie  following  symbols  are  used  in  this  report: 


A,B,C 

c* 

F 

f 

g 

g.l.b. 


constants 
critical  velocity 

perturbation  potential  in  xT]-plane 
function  of  x,t] 

thickness  distribution  along  chord 
greatest  lo-wer  bound 


Kc 


L  ,M,N,0 


M 

Mc 

^0 


n 

P 

CL 


s,t 


U 


functional  determinant ,  or  Jacobian 

transonic  similarity  parameter 

critical  value  of  similarity  parameter 

potential  limit  of  similarity  parameter 

constants,  identified  by  superscripts  and  subscripts 

local  Mach  number 

critical  Mach  number 

potential  limit  Mach  number 

free- stream  Mach  number 

number  of  approximations 

particular  integral 

function  of  s 

elliptic  coordinates 

free-stream  velocity 

cartesian  coordinates 


2275 
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* 


r 

A 

5 

T) 

T 


<P 

u>,p,cr 

Subscripts : 


Mo2(i  +  IA  Mq2) 

ratio  of  specific  heats 

Laplacian 

constant 

transformed  y- coordinate 
thickness  ratio 
perturbation  potential 
constants,  identified  by  subscripts 


x,y,Ti,s,t  indicate  partial  derivatives 

h, j,k,m,p,r,v  summation  indices 
Superscripts : 


X,n,v 


number  of  approximation 
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TABLE  III  -  THIRD  APPROXIMATION  FOR  ELLIPSE 
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TABLE  IT  -  FOURTH  APPROXIMATION  FOR  ELLIPSE 


(a)  A 

hj 


2 

3 

4 

5 

6 

— 

7 

8 

9 

10 

B 

0.14320 

0.09286 

-0.16301 

0.04965 

0.01268 

0.00436 

0.00179 

— 

0.00007 

a 

-0.21480 

-0.43335 

0.31663 

0.02354 

-0.02353 

-0.02820 

-0.00936 

— 

-0.00007 

a 

0.07160 

0.57662 

0.19904 

-0.52940 

0.06649 

0.02305 

0.03632 

0.00024 

_ 

-0.26551 

-0.75659 

0.80076 

-0.07488 

0.05956 

-0.07065 

— 

-0.00132 

_ 

0.02977 

0.49494 

-0.38547 

0.02823 

-0.10221 

0.06084 

— 

0.00262 

_ 

_ 

-0.09795 

0.04100 

— 

0.04808 

-0.01895 

— 

-0.00222 

a 

— 

0.00693 

— 

-0.00464 

— 

0.00069 

<»  *  hj 

\J 

h\ 

4 

1 

0.16250 

3 

-0.48749 

S 

0.48749 

7 

-0.16250 

<c)  '  hj 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

13 

0.22498 

-0.09547 

-0.02332 

0.01018 

-0.00018 

-0.00006 

— 

-0.00000 

B 

_ 

-0.08592 

-0.04018 

0.06871 

-0.00294 

0.00241 

0.00141 

0.00034 

— 

_ 

0.00682 

0.07203 

-0.05206 

0.01638 

-0.01209 

-0.00192 

-0.00186 

— 

-0.00001 

— 

■ 

BH 

— 

— - 

-0.01328 

-0.04347 

0.00673 

-0.01152 

-0.01829 

0.00942 

— 

0.00005 

— 

— 

— - 

— 

— 

0.00083 

0.01523 

-0.01377 

0.00125 

-0.00639 

0.00716 

-0.00132 

-0.00028 

— - 

— 

— 

— 

-- 

— 

-0.00187 

0.00085 

0.00134 

-0.00066 

— 

-0.00021 

0.00012 

— 

— 

— 

— 

— 

0.00009 

— 

-0.00008 

— 

0.00002 

— 

-0.00000 

(4)  M  h<} 
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to 

03 
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I 

r- 

1 

1 

1 

1 

1 

1 

r^- 

03 

o 

00 

1 

to 

in 

in 

CO 

rtf 

r- 

CO 

in 

LO 

03 

gjgjTO 

co 
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in 

rtf 

in 

i H 

o 

03 

O 

CO 

KISH 

o 

rH 

03 

03 

1 

o 

■ 

rH 

03 

CO 

r- 

o 

o 

in 

1 — 1 

rtf 

rtf 

to 

CO 

03 

to 

to 

CO 

in 

CO 

03 

rH 

rH 

rH 

rH 

o 

m 

CO 

rH 

f — 1 

o 

o 

H 

1 

m 

B 

03 

03 

■M 

m 

o 

O 

to 

o 

o 

to 

o 

o 

to 

CO 

Hgfl 

• 

o 

rH 

rH 

1 

I 

rH 

to 

in 

r- 

03 

rH 

rH 

TABLE  V  -  FIFTH  APPROXIMATION  FOR  ELLIPSE  -  Concluded 
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11 

1 

1 

1 

i 

i 

i 

i 

i 

i 

1 

| 

i 

i 

-0.00036 

10 

1 

1 

1 

i 

i 

i 

i 

i 

1 

t 

1 

i 

i 

i 

i 

i 

i 

o> 

1 

1 

1 

i 

t 

i 

t 

i 

i 

1 

1 

1 

0.01210 

i 

t 

i 

00 

CO 

o 

o 

o 

o 

o 

-0.00073 

LO 

* 

o 

o 

-0.02274 

98910*0“ 

0.00154 

r- 

i 

i 

i 

i 

1 

1 

-0.02476 

i 

i 

i 

i 

i 

CO 

0.00564 

to 

(—1 

CD 

8 

o 

l 

0.20414 

0.15737 

i 

i 

i 

i 

i 

i 

10 

0.00598 

0.00425 

-0.10345 

-0.09638 

0.01988 

-0.00123 

-0.02889 

0.18572 

0.14445 

-0.01313 

i 

i 

t 

i 

i 

to 

1 

2 

-0.14818 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

/ 

/a 

H 

to 

LO 

r- 

CD 

11 

* 
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M 
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TABLE  VI  -  SIXTH  APPROXIMATION  FOR  ELLIPSE 


(a)  A  hJ 


X 

2 

3 

4 

5 

6 

— 

7 

8 

9 

1 

-0.35512 

-0.20367 

-0.54907 

0.53760 

0.52224 

0.39212 

0.01983 

-0.00962 

3 

0.53268 

0.95047 

2.16610 

-1.13329 

-3.53878 

2.24494 

0.27821 

-0.05669 

5 

-0.17756 

-1.26381 

-3.67930 

-0.80210 

8.32920 

-3.82950 

-0.96660 

0.19693 

7 

0.58229 

3.23792 

4.50042 

-8.29361 

0.43150 

1.86413 

0.21831 

9 

— 

-0.06528 

-1.43991 

-5.29265 

1.98519 

6.03432 

-2.58489 

-0.95363 

11 

- _ 

... 

0.28437 

2.85636 

2.82542 

-7.73804 

2.34101 

0.50963 

13 

— 

— - 

-0.02010 

-0.75951 

-2.64252 

4.22655 

-1.11805 

0.50443 

15 

— 

_ 

— 

0.09828 

0.97957 

-1.11033 

0.22088 

-0.54201 

17 

. —  - 

-0.00511 

-0.18374 

0.13964 

-0.01486 

0.14955 

19 

— 

— 

— 

0.0X776 

-0.00697 

— 

-0.01773 

21 

— 

—  —md 

— 

— 

-0.00072 

— 

— 

0.00082 

6 

(a)  A  -  Concluded 


X  j 

h\ 

10 

11 

12 

13 

14 

15 

16 

18 

1 

0.00467 

-0.00511 

0.00088 

-0.00052 

-0.00014 

-0.Q0002 

-0.00000 

-0.00000 

3 

0.00801 

0.01975 

0.00047 

0.00045 

0.00043 

0.00006 

0.00002 

0.00000 

5 

-0.02668 

-0.05440 

-0.01173 

0.00009 

-0.00080 

-0.00001 

-0.00007 

-0.00000 

7 

-0.33773 

0.17216 

0.01806 

0.00866 

0.00131 

-0.00007 

0.00008 

0.00000 

9 

0.69461 

-0.17444 

0.00844 

-0.03344 

-0.00290 

-0.00018 

-0.00005 

-0.00000 

11 

0.11262 

-0.17842 

-0.02885 

0.03281 

0.00916 

0.00084 

0.00007 

0.00000 

13 

-1.00872 

0.39807 

-0.01228 

0.02266 

-0.01671 

-0.00085 

-0.00029 

0.00001 

15 

0.66684 

-0.19573 

0.05158 

-0.05514 

0.01389 

-0.00039 

0.00059 

-0.00003 

17 

-0.12088 

0.01812 

-0.03191 

0.02680 

-0.00423 

0.00115 

-0.00053 

0.00004 

19 

0.00725 

_  _  _ 

0.00568 

-0.00239 

_ _ 

-0,00057 

0.00017 

-0.00003 

21 

-0.00033 

_ — 

_ 

0.00005 

... 

0.00001 
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TABLE  VI  -  SIXTH  APPROXIMATION  FOR  ELLIPSE  -  Continued 


(*>  ®  hi 


\i 

h\ 

4 

5 

6 

7 

8 

9 

10 

12 

1 

-0.44175 

0.27651 

0.17244 

-0.21716 

-0.07427 

0.00325 

-0.01077 

-0.00071 

3 

1.32524 

0.07448 

-0.95724 

1.36008 

0.22528 

0.09495 

0.02864 

0.00142 

5 

-1.32524 

-2.37617 

1.17188 

-1.02543 

-0.60779 

-0.40598 

-0.11295 

-0.00071 

7 

0.44175 

3.51929 

1.55575 

-3.43737 

1.21433 

0.34626 

0.45626 

0.00399 

9 

-1.66999 

-3.93061 

5.82295 

-1.13510 

0.47585 

-0.84175 

-0.02074 

11 

0.17588 

2.42446 

-2.77921 

0.37754 

-0.89299 

0.69364 

0.02897 

13 

-0.46759 

0.27615 

— 

0.41725 

-0.20907 

-0.03169 

15 

— 

— 

0.03092 

— 

— 

-0.03859 

— 

0.00947 

\J 

h.\ 

1 

-0.38239 

3 

1.24277 

5 

-1.43397 

7 

0.66918 

9 

-0.09560 

cc*  to 
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TABLE  VI  -  SIXTH  APPROXIMATION  FOR  ELLIPSE  -  Continued 

(a)  l  hJ 


V  j 

h\ 

1 

2 

3 

4 

5 

6  ’ 

7 

8 

9 

10 

1 

0.34464 

-0.11837 

-0.02546 

-0.05231 

0.02720 

0.01467 

-0.00949 

-0.00061 

-0.00011 

0.00003 

3 

_ 

-0.10654 

-0.40107 

0.52581 

-0.06792 

-0.13205 

0.06802 

0.00625 

-0.00046 

0.00016 

5 

_ 

0.00846 

0.07899 

0.27792 

-0.66704 

0,58684 

-0.18449 

-0.03118 

0.00119 

-0.00076 

7 

_ 

_ 

-0.01456 

-0.09487 

-0.12642 

0.65280 

-0.67976 

0.21063 

0.01291 

-0.00314 

9 

_ 

_ _ 

0.00091 

0.02215 

0.08919. 

-0.06701 

-0.10896 

0.08921 

0.00838 

-0.01008 

11 

_ 

_ 

' 

-0.00271 

-0.02956 

-0.02921 

0.09997 

-0.03921 

-0.02279 

0.02609 

13 

_ 

_ 

_ 

0.00013 

0.00527 

0.01955 

-0.03495 

0.01065 

-0.00476 

0.01374 

15 

... 

... 

-0.00049 

-0.00517 

0.00631 

-0.00137 

0.00373 

-0.00533 

17 

— — 

0,00002 

0.00073 

-0.00058 

0.00007 

-0.00064 

0.00063 

19 

... 

-0.00005 

0.00002 

0.00006 

-0.00005 

21 

_ 

_ 

-__r 

_ 

0.00000 

_ 

_ 

-0.00000 

6 

(d)  L  ^  -  Concluded 


h\ 

11 

12 

13 

14 

15 

16 

17 

18 

19 

21 

1 

-0.00004 

0.00001 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

_ 

-0.00000 

... 

3 

0.00018 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

5 

-0.00057 

-0.00010 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

7 

0.00239 

0.00020 

0.00007 

0.00001 

-0.00000 

0.00000 

_ _ 

0.00000 

... 

9 

-0,00436 

0.00001 

-0.00038 

-0.00003 

-0.00000 

-0.00000 

_ 

-0.00000 

11 

-0.01128 

0.00051 

0.00068 

0.00012 

0.00001 

0.00000 

0.00000 

13 

-0.00829 

-0.00073 

0.00227 

-0.00062 

-0.00002 

-0.00000 

- 

0.00000 

15 

0.00188 

-0.00060 

0.00098 

-0.00048 

0.00001 

0.00002 

-0.00000 

_ _ 

_ 

17 

-0.00011 

0.00022 

-0.00022 

0.00005 

-0.00002 

0.00002 

-0.00000 

0.00000 

_ 

19 

_ 

-0.00003 

0.00001 

_ 

0.00000 

-0.00000 

_ 

0.00000 

0.00001 

_ 

21 

-i —  - 

0.00000 

— 

— 

-0.00000 

— 

-0.00000 

— 

0.00000 

TABLE  VI  -  SIXTH  APPROXIMATION  FOR  ELLIPSE  -  Concluded 
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-0.00417 
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TABTE  IX  -  THIRD  APPROXIMATION  FOR  KAPLAN  SECTION 
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-0.10413 


TABLE  X  -  FOURTH  APPROXIMATION  FOR  KAPLAN  SECTION 


4 


\  J 

2 

3 

4 

5 

6 

7 

8 

9 

1 

-0.27369 

-0.10158 

-0.70164 

0.16384 

-0.25999 

-0.08941 

0.27562 

0.00827 

3 

0.41053 

0.47402 

2.21465 

0.20698 

-0.14470 

0.34639 

-0.94946 

0.74600 

5 

-0.13684 

-0.63029 

-2.54101 

-2.37393 

-0.09017 

-0.37686 

3.50541 

-2.30284 

7 

— 

0.29040 

1.28382 

3.50709 

2.64389 

-6.52365 

1.77718 

-2.32167 

9 

--- 

-0.03256 

-0.31320 

-1.68637 

-3.74815 

-3.76360 

2.16703 

1.62381 

11 

— 

—  - 

0.06195 

0.18240 

1.95475 

2.08185 

4.82315 

-3.15023 

13 

--- 

— 

-0.00437 

— 

-0.38139 

-0.24116 

-4.05218 

1.40343 

15 

— 

— — — 

— 

— 

0.02576 

--- 

0.80760 

-0.10562 

17 

— 

— 

— 

— 

--- 

--- 

-0.05443 

— 

19 

— 

— 

— 

— 

— 

--- 

— 

— 

21 

--- 

— 

--- 

— - 

— 

— 

. . . . .  j 

— - 

4 

(a)  A  ^  -  Concluded 


m 

10 

11 

12 

13 

14 

15 

16 

18 

1 

0.30919 

-0.13252 

-0.03337 

-0.00319 

-0.03780 

0.00311 

-0.00705 

-0.00035 

-1.56053 

0.63407 

-0.17619 

0.30166 

0.04789 

0.00980 

0.01409 

0.00035 

5 

1.92675 

-0,67733 

0.58936 

-0.64026 

-0.13063 

-0.06476 

-0.00705 

_ 

7 

1.37105 

1.37105 

-0.18288 

0.12860 

0.42178 

0.09728 

-0.02032 

_ 

m 

5.14985 

-2.08886 

-0.71016 

-0.22855 

-0.46413 

0.00347 

0.07111 

-0.00280 

11 

-2.88950 

-2.33567 

0.76885 

1,04169 

0.41837 

-0.03373 

-0.07825 

0.00796 

■ 

-0.99893 

0.85868 

-0.95828 

-0.07389 

0.02888 

-0.00818 

m 

1.55071 

-2.77829 

1.12471 

-0.09933 

-0.01642 

0.00786 

17 

-0.71285 

0.22345 

-1.01520 

-0.42191 

0.32418 

0.02302 

-0.01576 

19 

0.04935 

— 

0.22411 

—  - 

-0.18607 

-0.00802 

0.01698 

21 

— 

— 

-0.01630 

— 

— - 

0.01994 

— 

-0.00606 

4 

00  B 


X 

4 

6 

8 

10 

12 

1 

0.01714 

0.02471 

0.02657 

0.00976 

— 

3 

-0.05142 

-0.11585 

-0.07267 

-0.01881 

-0.00395 

5 

0.05142 

0.24496 

0.17382 

0.07426 

0.01578 

7 

-0.01714 

-0.22984 

-0.35115 

-0.09793 

-0.01973 

9 

— 

0.07599 

0.33867 

0.18995 

0.01776 

11 

— 

— 

-0.11524 

-0.19764 

-0.03354 

13 

— 

— 

0.07041 

0.03749 

15 

— 

— 

— 

-0.01381 
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TABLE  X  -  FOURTH  APPROXIMATION  FOR  KAPLAN  SECTION  -  Concluded 


4 

(c)  L 


\J 
h \ 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

0.37947 

-0.09123 

-0.01270 

-0.04618 

0.00683 

-0.00719 

-0.00186 

0.00450 

0.00010 

0.00314 

3 

--- 

-0.08211 

-0.36592 

0.30799 

0.01294 

-0.00727 

0.00866 

-0.01764 

0.01036 

-0.01727 

5 

_ 

0.00652 

0.03939 

6.28741 

-0.30613 

0.01610 

-0.01570 

0.09171 

-0.04112 

0.02595 

7 

_ 

-0.00726 

-0.03903 

-0.14613 

-0.21970 

-0.78172 

-0.45988 

0.05554 

-0.04627 

9 

_ _ 

0.00045 

0.00482 

0.03011 

0.08374 

0.11761 

-0.10872 

-0.30580 

0.27104 

11 

... 

_ 

-0.00059 

-0.00190 

-0.02300 

-0.02891 

-0.08519 

0.07876 

0.12864 

13 

_ 

_ 

0.00003 

— 

0.00287 

0.00201 

0.03859 

-0.01595 

0.04919 

15 

... 

_ 

_ 

— 

— 

-0.00014 

— 

-0.00502 

0.00073 

-0.02738 

17 

_ 

_ 

— 

... 

— 

0.00024 

... 

0.00377 

19 

_ 

_ _ 

... 

— 

... 

... 

-0.00019 

21 

... 

... 

— 

— 

... 

— 

... 

... 

... 

... 

(c)  L  hj  -  Concluded 


h\ 

11 

12 

13 

14 

15 

16 

17 

18 

19 

21 

1 

-0.00110 

-0.00023 

-0.00002 

-0.00019 

0.00001 

-0.00003 

.... 

-0.00000 

— 

— 

3 

0.00566 

-0.00130 

0.00189 

.0,00026 

0.00005 

0.00006 

— 

0.00000 

— 

— 

5 

-0.00706 

0.00498 

-0.00445 

-0.00076 

-0.00082 

-0.00011 

— 

— 

— 

— 

7 

0.01904 

-0.00198 

0.00107 

0.00287 

0.00035 

-0.00010 

— 

— 

— 

— 

9 

-0.05222 

-0.01127 

-0.00260 

-0.00404 

0.00002 

0.00041 

— 

-0.00001 

— 

— 

11 

-0.14319 

0.03191 

0.02170 

0.00558 

-0.00032 

-0.00058 

0.00004 

— 

13 

-0.12831 

0.04140 

0.05407 

-0.03549 

-0.00132 

0.00033 

— 

-0.00005 

— 

— 

15 

0.03032 

-0.01919 

0.04961 

-0.03878 

0.00756 

-0.00053 

’  - 

0.00008 

— 

— 

17 

-0.00133 

0.00700 

-0.01122 

0.00454 

-0.00507 

-0.00070 

0.00404 

-0.00045 

— 

— 

19 

-0.00103 

0.00014 

0.00137 

0.00008 

— 

-0.00046 

-0.00006 

— 

21 

... 

0.00008 

— 

— 

-0.00009 

— 

— 

0.00005 

— 

-0.00002 

(4)  M  hJ 


3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

15 

1 

0.00114 

_ 

0.00071 

— 

0.00042 

— 

0.00010 

— 

— 

— 

— 

3 

-0.07900 

-0.00735 

-0.00429 

_ 

-0.00132 

— 

-0.00054 

— 

-0.00003 

— 

— 

5 

-0.00571 

0.23739 

0.02227 

_ 

0.00446 

— 

0.00099 

— 

0.00013 

— 

— 

7 

0.00052 

_ _ 

0.01768 

-0.14591 

-0.02341 

— 

-0.00192 

— 

-0.00021 

— 

— 

g 

-0.00169 

_ 

-0.01992 

-0.09021 

0.01000 

— 

0.00028 

— 

— 

11 

_ 

_ 

0.00202 

0.00941 

0 . 10617 

-0.00146 

— 

— 

13 

_ 

_ 

_ 

— 

— 

— 

-0.00102 

— 

-0.00150 

-0.03303 

— 

15 

— 

_ — 

— 

— 

— 

— 

— 

— 

0.00017 

_ 

0.00331 
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Figure  1.  -  Reciprocal  of  transf ormation  coefficient. 


Perturbation  velocities ,  <pv  and 
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Figure  2.  -  Parabolic  curve  of  transonic  equation. 


42 


MCA  TIT  2511 


Figure  4.  -  Functional  determinant  of  transformation 
x  =  cosh  s  cos  t,  r\  =  sinh  s  sin  t. 


Figure  5.  -  Kaplan  section  and  region  of  validity  of  simplified  flow  equation.  Thickness  ratio,  0.1;  free-stream 

Mach  number,  0.9. 
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-stream  Mach  number 


Ratio  of  successive  terms 
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Number  of  approximations,  n 


Figure  8.  -  Estimation  of  potential  limit  similarity  parameter  for 

elliptic  cylinder. 
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x  =  cos  t 


Figure  13.  -  Telocity  increment  on  surface  of  Kaplan  section.  Thickness 
ratio,  0.10;  free -stream  Mach  number,  0.75;  three  approximations. 
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Number  of  approximations/  n 


Figure  14.  -  Estimation  of  potential  limit, 
similarity  parameter  for  Kaplan  section. 
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Figure  15.  -  Maximum  Telocity  increment  on  surface  of  Kaplan  section. 


